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Abstract 

The Configuration Interaction method is applied to investigate the possibility of 
positron binding to the metastable beryllium {ls 2 2s2p 3 P°) state. The largest cal- 
culation obtained an estimated energy that was unstable by 0.00014 Hartree with 
respect to the Ps + Be + (2s) lowest dissociation channel. It is likely that positron 
binding to parent states with non-zero angular momentum is inhibited by centrifugal 
barriers. 
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1 Introduction 



In the last several years there has been substantial progress studying the 
physics of positron binding to atoms. As late as 1997, there was no rigorous 
evidence for existence of positron-atom bound states [1], although the pos- 
sibility of positron binding had been often invoked to explain the existence 
of very large annihilation rates for positrons annihilating in gases [2,3,4,5,6]. 
Towards the end of 1997, positron binding to atomic lithium was established 
in two independent calculations [7,8]. Following these initial studies, positron 
binding to a further 10 elements of the periodic table was demonstrated [9]. 
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One of the features of these positronic atoms is that all the atoms binding a 
positron have ionization potentials close to 6.80 eV (the Ps binding energy) 
and the binding energies are generally larger for the atoms with their ioniza- 
tion energies closest to 6.80 eV. 

A related question is whether positrons can form bound states with excited 
electronic states. There is one excited atomic state known to bind a positron, 
namely the metastable He( 3 S e ) state [10,11]. The stability of e + He( 3 S e ) state 
was established with the Stochastic Variational Method (SVM) [12,13,7] fol- 
lowing an initial inconclusive study by Drachman et al [14]. 

The answer to this question is relevant to positron- atom collision physics [15], 
since such states may manifest themselves as Feshbach resonances in the elastic 
or excitation cross sections. There has been a suggestion that the rich resonant 
structures prominent in electron- atom scattering [16] are largely absent from 
the positron-atom spectrum [17,9,15]. 

This work investigates the possibility of positron binding to the low-lying 
triplet odd state of neutral beryllium, i.e. the Be(ls 2 2s2p 3 P°) system. This 
state has been identified as one of the most likely excited states to bind a 
positron [18,9]. This conjecture has largely been based upon the fact that 
its ionization potential (IP) of 6.598 eV [19] is quite close to the Ps binding 
energy. This IP is between that of magnesium and calcium, both of which 
bind a positron with a binding energy of 0.4 eV [9,20,21]. Furthermore, 
the dipole polarisability of Be( 3 P°) is a d = 39.02 a 3 , [22], greater than that of 
ground state Be (a d = 37.7 a 3 , [23]). 

This metastable state warrants investigation for another reason; the parent 
atom has an angular momentum greater than zero with the loosely bound 2p 
electron carrying the angular momentum. The quantitative impact that the 
angular momenta has on the ability to bind a positron is largely unknown. 
Angular momentum coupling considerations suggests that the existence of a 
repulsive centrifugal barrier will act to inhibit positron binding by making it 
less energetically favorable to form a Ps-cluster. Either, the electron-positron 
pair will be in a relative p-state with a smaller energy gain than the Ps(ls) 
state, or the total angular momentum of the Ps cluster relative to the nucleus 
will be greater than zero. 



2 Details of the Calculation 

Since the ionization energy of the Be(2s2p 3 P°) system is just less than 6.80 eV, 
the threshold for a stable positron complex is that for the Be + (2s) + Ps(ls) 
dissociation channel. This energy was —0.9192086 Hartree for the present 
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model potential. An energy level diagram showing the ground and low-lying 
metastable states of Be, and some negative ion and positronic ion energies are 
shown in Figure 1. 

The CI method as applied to positron-atomic systems with two valence elec- 
trons and a positron has been discussed previously [24,25,20], but a short de- 
scription is worthwhile. The model Hamiltonian is initially based on a Hartree- 
Fock (HF) wave function for the neutral atom ground state. One- and two- 
body semi-empirical polarization potentials are added to the potential field of 
the HF core and the parameters of the core-polarization potentials defined by 
reference to the spectrum of Be + [26]. 

All calculations were done in the frozen-core approximation. The general ef- 
fective Hamiltonian for the system with N e valence electrons and a positron 

was 



1 JV e 1 N e 1 N e 

H = -^l ~ E + E — + E(^(r.) + V exc ( ri ) + V pl (Ti)) 
z i=\ z i<j r *j 1=1 

- E — - ^Kr(ro) + V pl (r ) - E V p^ r i) + E V P^ii r o) ■ (1) 

The direct potential (Vdi r ) represents the interaction with the HF Is 2 electron 
core. The direct part of the core potential is attractive for electrons and re- 
pulsive for the positron. The exchange potential (V exc ) between the valence 
electrons and the HF core was computed without approximation. 

The one-body polarization potential (V p i) was a semi-empirical polarization 
potential derived from an analysis of the Be + spectrum. It has the functional 
form 

<*d9 2 (r) 
2r 4 

The factor a a is the static dipole polarizability of the core and g 2 (r) is a cutoff 
function designed to make the polarization potential finite at the origin. The 
same cutoff function has been adopted for both the positron and electrons. In 
this work, g 2 (r) was defined to be 



VAr) = -^kr 1 - (2) 



g 2 (r) = 1 - exp(-r 6 /p 6 ) , (3) 

where p is an adjustable cutoff parameter. The core dipole polarizability was 
set to 0.0523 Oq while p was set to 0.95 a [25]. The two-body polarization 
potential (V p2 ) is defined as 

V p2 (ri, Tj) = -3-3 (r< • r j )^(r i ) 5 (r j ) . (4) 

i 3 
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This model potential gives a Be(2s2p 3 P°) binding energy of 0.91147 Hartree 
with respect to the Be 2+ threshold. The experimental binding energy with 
respect to this threshold is 0.91888 Hartree [27]. 

The CI basis was constructed by letting the two electrons and the positron 
form all the possible total angular momentum Lt = 1 configurations, with the 
two electrons in a spin-triplet state, subject to the selection rules, namely 



In these rules £q, £\ and £2 are respectively the orbital angular momenta of 
the positron and the two electrons. 

Our two-electron-positron calculations with non-zero total angular momen- 
tum were first validated against the previous Lj- = 1 and 2 PsH calculations 
of Tachikawa [28] . Using their Gaussian-type orbitals we reproduced their re- 
ported energy and annihilation rates. Note that the PsH states with L T — 1 
and 2 are unbound [29,28,30]. 

For the e + Be( 3 P°) calculations, the Hamiltonian was diagonalized in a CI basis 
constructed from a very large number of single particle orbitals, including 
orbitals up to £ = 12. There was a minimum of 12 radial basis functions for 
each £. The largest calculation was performed with L max = 12 and L int = 3 and 
gave a CI basis dimension of 498750. The resulting Hamiltonian matrix was 
diagonalized with the Davidson algorithm [31], and a total of 1376 iterations 
were required for the largest calculation. 



3 Results 

Various e + Be( 3 P°) expectation values, such as energy and mean distance of 
the electron and positron from the nucleus are given in Table 1. The 27 annihi- 
lation rate [9] for annihilation with the core and valence electrons are denoted 
r c and r„ respectively. The calculations shown in Table 1 have a minimum of 
12 radial basis functions for each £ with L max = 12 and L int = 3. The largest 
calculation remains unbound by —0.0021904 Hartree. 

To demonstrate that L int = 3 is sufficient for this system, also shown in Table 
1 is the result of a L max = 12 and L mt = 2 calculation. The L max parameter 
needs to be large since it determines the extent to which electron-positron 
correlations are incorporated into the wavefunction. However, L int is largely 
concerned with electron-electron correlations and, for example, setting L int = 3 



max(£ ,£ 1 ,£ 2 ) < 
min(£i,£ 2 ) < 



(5) 
(6) 
(7) 
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for the PsH and e + Be ground states recovered 99.4% and 98.5% respectively 
of the Ps and positron binding energies for a given L max [25]. 

The main problem afflicting CI calculations of positron-atom interactions is 
the slow convergence of the expectation values with L max [32,33,9,34]. One 
way to determine the L max — > oo expectation values is to assume that the 
successive increments, AXl, to any expectation value (X) scale as 1/L P as 
L increases [32]. However, arguments based on 2nd-order perturbation theory 
suggests that the energy increments scale asymptotically as 1/(1/ + \) PE ~ 
1/(L + |) 4 [35,36]. While existing CI calculations are consistent with this idea 
[25,37,34,38], the actual size of the exponent, pe is significantly smaller than 
4 at L max pa 10. This can be seen in Table 1, where the power law factor pe 
relating successive i?(L max — 2), E(L max — l), and E(L max ) are tabulated. Ap- 
plying a simple power law extrapolation with pe = 2.078 leads to a prediction 
of binding (E = —0.9207021), but not much credence should be placed in this 
since fixing pe at its L max = 12 value of 2.078 will result in an extrapolation 
that will overestimate the contribution of the orbitals with L > L max [39,34]. 

An improved method, borrowing from existing work in atomic structure physics 
[40,41,21], is to estimate the L max — > oo limit for the energy by assuming that 

. „ Ae Be Ce / x 

The factors Ae, Be and Ce are determined from four calculations at suc- 
cessively larger values of L max . Applying eq. (8) to the data in Table 1 at 
L max = 12 one finds A E = -34.601, B E = 477.71 and C E = -1829.3, with 
E = —0.919070 Hartree as our present best energy estimate. This is 0.0001390 
Hartree higher the threshold for binding at —0.9192086 Hartree. 

Usage of an inverse power series restricted to the first two terms of eq. (8) 
results in A E = -21.875 and B E = 172.29, and an e+Be( 3 P°) energy of 
E = —0.9188267, which is also unbound. As has been seen in other CI calcu- 
lations of positron-atom systems, even when the three-term and the two-term 
coefficents are significantly different, their extrapolated energies can still lie 
relatively close together 



There are two aspects where convergence can be incomplete. Besides the size of 
L max , the number of LTOs for each value of I still needs to be further increased. 
This is best seen in Figure 2 which depicts the energies as a function of L max . 
Two sets of e + Be( 3 P°) data are shown, the first shows the e + Be( 3 P°) energies 
from Table 1 which used 12 LTOs, and their extrapolated estimates using 
eq. (8). The second set shows the results from initial exploratory calculations 
also with L max = 12 and L int = 3, but only with a minimum of 8 LTOs 
per £. Even a cursory glance at Figure 2 suggests that the CI prognosis for 
variationally establishing the existence of a bound state is not promising. Also 
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shown in Figure 2 are the Be( 3 P°) energies using the underlying e + Be( 3 P°) 
basis without the positron (the 12 LTOs basis). 

Even though the valence annihilation rate from the largest calculation is only 
0.348 x 10 9 sec -1 , the slow convergence of T v with L max suggests the presence 
of a well defined Ps cluster. Fitting the four largest calculations in Table 1 
to the three-term form 

results in an extrapolated annihilation rate of T v = 0.805 x 10 9 sec -1 . However, 
this value should be regarded as only notional. After all, if this system is indeed 
unbound, r„ should approach Tp s = 2.008 x 10 9 sec -1 in the limit of an infinite 
radial basis [42]. 

It is also noticed that (r p ) in Table 1 steadily decreases as L max increases. This 
also should not be taken as an indicator of binding since the behavior of (r p ) 
with L max is not straightforward. For example, the PsH (r p ) decreases when 
L max changes from to 3, but then starts increasing as L max increases from 3 to 
9 [25,34]. Indeed, this is a trend noticed in CI calculations of positronic atoms 
with IP < 6.8eV (eg. e+Li [43], e+Ca and e+Sr [20,38]), which have shown an 
initial (r p ) decrease to a minimum, which then increases as L max — > oo. The 
behavior of (r p ) with L max for an unbound system is unknown. 



4 Conclusions 



A large-scale CI calculation of the e + Be( 3 P°) ground state has been performed. 
The largest calculation gave an energy that was 0.00219 Hartree above the 
threshold for binding. Using an extrapolation method to estimate the contri- 
bution from the higher-partial waves gave an energy that was 0.00014 Hartree 
above threshold. By its very nature the present calculation is unable to give 
definitive proof that the e + Be( 3 P°) system does not have a bound state. While 
a variational calculation of the present kind can be used to give proof of bind- 
ing, it cannot be used to give proof of non-binding. A converged calculation 
of low energy Be + -Ps(ls) scattering in the 4 P° channel would be needed to 
establish lack of binding. 

What can be inferred from the present calculation is that a e + Be( 3 P°) bound 
state (assuming one existed) would have a structure similar to positronic 
lithium, i.e. it would consist of a Ps-like object weakly bound to the Be + (2s) 
core. The size of the CI expansion required to establish binding would be enor- 
mous. For example, it is necessary to include orbitals with L max 30 to get 
an energy that is lower than the Ps + Li + threshold [43]. Another conclusion 
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is that the existence of centrifugal barriers does seem to inhibit the binding 
of positrons to parent atoms with non-zero angular momentum. 



Acknowledgments 

The authors would like to thank Shane Caple of CDU for computer support. 
Some preliminary calculations were performed on computing facilities made 
possible by the Research Corporation. One author (MB) would like to thank 
Dr. Brett Esry for discussions on three-body systems with non-zero angular 
momentum. The authors would also like to thank Dr. Masanori Tachikawa 
for access to unpublished data on positronic systems with non-zero angular 
momentum. 



References 



[1] D. M. Schrader, Bound States of Positrons with Atoms and Molecules: Theory, 
Nucl. Instrum. Methods Phys. Res. B 143 (1998) 209. 

[2] D. A. L. Paul, L. Saint-Pierre, Rapid Annihilations of Positrons in Polyatomic 
Gases, Phys. Rev. Lett. 11 (1963) 493. 

[3] H. C. Khare, P. R. Wallace, G. G. Bach, A. Chodos, Annihilation of positrons 
in liquid helium and the positron-helium bound state, Can. J. Phys. 42 (1964) 
1522. 

[4] V. I. Goldanskii, Y. S. Sayasov, On the resonance annihilation of positrons in 
collisions of neutral atoms or molecules, Phys. Lett. 13 (1964) 300. 

[5] J. D. McNutt, D. A. Johnson, V. B. Summerour, Annihilations of Positrons in 
Methane, Phys. Rev. B 4 (1971) 36. 

[6] C. M. Surko, A. Passner, M. Leventhal, F. J. Wysocki, Bound states of positrons 
and large molecules, Phys. Rev. Lett. 61 (1988) 1831. 

[7] G. G. Ryzhikh, J. Mitroy, Positronic lithium, an electronically stable Li-e + 
ground state, Phys. Rev. Lett. 79 (1997) 4124. 

[8] K. Strasburger, H. Chojnacki, Quantum chemical study of simple positronic 
systems using explicitly correlated Gaussian functions - PsH and PsLi + , 
J. Chem. Phys. 108 (1998) 3218. 

[9] J. Mitroy, M. W. J. Bromley, G. G. Ryzhikh, Positron and positronium binding 
to atoms and ions, J. Phys. B 35 (2002) R81. 

[10] G. G. Ryzhikh, J. Mitroy, A metastable state of positronic helium, J. Phys. B 
31 (1998) 3465. 



7 



[11] J. Mitroy, Expectation values of the e + He( 3 S e ) system, Phys. Rev. A 72 (2005) 
032503. 

[12] K. Varga, Y. Suzuki, Precise solution of few-body problems with the stochastic 
variational method on a correlated Gaussian basis, Phys. Rev. C 52 (1995) 2885. 

[13] Y. Suzuki, K. Varga, Approach to Quantum-Mechanical Few-Body Problems, 
172, Springer, New York, 1998. 

[14] R. J. Drachman, Y. K. Ho, S. K. Houston, Positron attachment to helium in 
the 3 5 state, J. Phys. B 9 (1976) L199. 

[15] C. M. Surko, G. F. Gribakin, S. J. Buckman, Low-energy positron interactions 
with atoms and molecules, J. Phys. B 38 (2005) R57. 

[16] S. J. Buckman, C. W. Clark, Atomic negative-ion resonances, Rev. Mod. Phys. 
66 (1994) 539. 

[17] J. P. Sullivan, S. J. Gilbert, S. J. Buckman, C. M. Surko, Searches for resonances 
in the scattering of low-energy positrons from atoms and molecules, J. Phys. B 
34 (2001) L467. 

[18] D. M. Schrader, Antimatter Compounds, in: C. M. Surko, F. A. Gianturco 
(Eds.), New Directions in Antimatter Physics and Chemistry, Kluwer Academic 
Publishers, The Netherlands, 2001, p. 263. 

[19] S. Bashkin, J. O. Stoner, Atomic Energy Levels and Grotrian Diagrams, North 
Holland, Amsterdam, 1975. 

[20] M. W. J. Bromley, J. Mitroy, Configuration Interaction calculations of positron 
binding to group II elements, Phys. Rev. A 65 (2002) 062505. 

[21] M. W. J. Bromley, J. Mitroy, Convergence of partial wave expansion of the He 
ground state (2006) in preparation. 

[22] J. Mitroy, M. W. J. Bromley, Properties of the triplet metastable states of 
alkaline-earth-metal atoms, Phys. Rev. A 70 (2004) 052503. 

[23] J. Mitroy, M. W. J. Bromley, Semi-empirical calculation of van der Waals 
coefficients for alkali-metal and alkaline-earth-metal atoms, Phys. Rev. A 68 
(2003) 052714. 

[24] M. W. J. Bromley, J. Mitroy, G. G. Ryzhikh, Configuration interaction 
calculations of positronic atoms and ions, Nucl. Instrum. Methods Phys. Res. B 
171 (2000) 47. 

[25] M. W. J. Bromley, J. Mitroy, Configuration Interaction calculations of PsH and 
e+Be, Phys. Rev. A 65 (2002) 012505. 

[26] D. W. Norcross, M. J. Seaton, Energy levels for Be I calculated using a model 
potential and cores approximations, J. Phys. B 9 (1976) 2983. 

[27] NIST Atomic Spectra Database Version 3.0 (2005). 

URL http : //physics . nist . gov/cgi-bin/AtData/main$_$asd 

8 



[28] M. Tachikawa, Simultaneous optimization of Gaussian type function exponents 
for electron and positron with full-CI wavefunction - application to ground and 
excited states of positronic compounds with multi-component molecular orbital 
approach, Chem. Phys. Lett. 350 (2001) 269. 

[29] D. Bressanini, M. Mella, G. Morosi, Stability and positron annihilation of 
positronium hydride L = 0, 1, 2 states: A quantum Monte Carlo study, 
Phys. Rev. A 57 (1998) 1678. 

[30] S. L. Saito, Multireference configuration interaction calculations of some low- 
lying states of positronium hydride, J. Chem. Phys. 118 (2003) 1714. 

[31] A. Stathopolous, C. Froese Fischer, A davidson program for finding a 
few selected extreme eigenpairs of a large, sparse, real, symmetric matrix, 
Comput. Phys. Commun. 79 (1994) 268. 

[32] J. Mitroy, G. G. Ryzhikh, Configuration interaction calculation of the Cue 4 " 
ground state, J. Phys. B 32 (1999) 2831. 

[33] V. A. Dzuba, V. V. Flambaum, G. F. Gribakin, C. Harabati, Calculation of the 
positron bound state with the copper atom, Phys. Rev. A 60 (1999) 3641. 

[34] J. Mitroy, M. W. J. Bromley, Convergence of CI single center calculations of 
positron-atom interactions (2006) in preparation. 

[35] C. Schwartz, Importance of Angular Correlations between Atomic Electrons, 
Phys. Rev. 126 (1962) 1015. 

[36] G. F. Gribakin, J. Ludlow, Convergence of partial-wave expansions for energies 
scattering amplitude and positron annihilation rates, J. Phys. B 35 (2002) 339. 

[37] M. W. J. Bromley, J. Mitroy, Variational calculation of positron-atom scattering 
using configuration-interaction-type wave functions, Phys. Rev. A 67 (2003) 
062709. 

[38] M. W. J. Bromley, J. Mitroy, Large dimension Configuration Interaction 
calculations of positron binding to the group II atoms, Phys. Rev. A 72 (2005) 
in prep. 

[39] M. W. J. Bromley, J. Mitroy, Configuration Interaction calculations of positron 
binding to Zn and Cd, Phys. Rev. A 65 (2002) 062506. 

[40] D. P. Carroll, H. J. Silverstone, R. P. Metzger, Piecewise polynomial 
configuration interaction natural orbital study of Is 2 helium, J. Chem. Phys. 
71 (1979) 4142. 

[41] R. N. Hill, Rates of convergence and error estimation formulas for the Rayleigh- 
Ritz variational method, J. Chem. Phys. 83 (1985) 1173. 

[42] J. Mitroy, M. W. J. Bromley, Comment on: Multi-reference configuration 
interaction calculations for positronium halides [J. Chem. Phys. 122, 054302 
(2005)], J. Chem. Phys. 123 (2005) 017101. 



9 



[43] M. W. J. Bromley, J. Mitroy, Positron and positronium interactions with Cu, 
Phys. Rev. A 66 (2002) 062504. 



[44] T. Andersen, H. H. Andersen, P. Balling, P. Kristensen, V. V. Petrunin, 
Structure and dynamics of the negative alkaline-earth ions, J. Phys. B 30 (1997) 
3317. 

[45] J. Mitroy, G. G. Ryzhikh, Improved description of LiPs, e + Be, NaPs and e + Mg, 
J. Phys. B 34 (2001) 2001. 

^ Be 2+ threshold 

-0.6 - 

Be + 2s( 2 S e ) 









Be 2p 


2 ( 3 P 6 ) Be-2p 3 ( 4 S°) 






6.598 eV 6.803 eV 




Be 2s2p( 3 P°) Be - 2s 2p 2 ( 4 P e ) Be + 


-Ps 


Be2s 2 ( 1 S e ) e + 


Be 



Fig. 1. The energy levels of electron and positron binding to neutral beryllium (in 
units of Hartree relative to the Be 2+ threshold). The experimental atomic binding 
energies are taken from the NIST compilation [27], the metastable negative ion 
experimental binding energies are from [44], and the theoretical e + Be binding energy 
is taken from a frozen-core SVM calculation [45]. The threshold for e + binding to 
Be( 3 P°) is denoted as Be++Ps. 
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Table 1. Results of CI calculations for e + Be( 3 P°) for a series of L max , with fixed L mt = 3 (the row denoted 12* was computed with 
Li n t = 2). The total number of electron and positron orbitals are identical and denoted by N. The 3-body energy of the e + -Be( 3 P°) 
system, relative to the energy of the Be 2+ core, is denoted by E (in Hartree). The threshold for binding is -0.9192086 Hartree, and e gives 
binding energy (in Hartree) against dissociation into Ps + Be + (2s). The mean electron-nucleus distance (r e ), the mean positron-nucleus 
distance (r p ), and the mean electron-positron separation (r| p ) are given in units of ao- The T v and T c columns give the valence and core 
annihilation rates respectively (in units of 10 9 sec -1 ). The pe column gives the power-series exponents from the L max , L max — 1, L max — 2 
energies. The results in the row oo used eq. (8) and eq. (9) to estimate the L max — ► oo correction. 
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Fig. 2. The energy (in units of Hartree) of e + -Be( 3 P°) (crosses) and extrapolated 
energies (diamonds) as a function of L max for two different calculations (with a mini- 
mum of 8 and 12 LTOs per i respectively). The Be( 3 P°) two-electron energy is shown 
by the squares for the 12 LTOs basis. The threshold for binding at E(Be + )+E(Ps) 
is shown as solid line. 
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